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Abstract 

The Henon-Heiles Hamiltonian was introduced in 1964 [jl|] as a mathematical model to describe 
the chaotic motion of stars in a galaxy. By canonically transforming the classical Hamiltonian to a 
Birkhoff-Gustavson normalform Delos and Swimm obtained a discrete quantum mechanical energy 
spectrum. The aim of the present work is to first quantize the classical Hamiltonian and to then diag- 
onalize it using different variants of flow equations, a method of continuous unitary transformations 
introduced by Wegner in 1994 [^]. The results of the diagonalization via flow equations are com- 
parable to those obtained by the classical transformation. In the case of commensurate frequencies 
the transformation turns out to be less lengthy. In addition, the dynamics of the quantum mechanical 
system are analyzed on the basis of the transformed observables. 

PACS-numbers: 03.65.-w (Quantum mechanics), 

05.45.H-b (Theory and models of chaotic systems) 

Keywords: Henon-Heiles Hamiltonian, Quantum chaos. Flow equations 

1 Introduction 

The Henon-Heiles Hamiltonian describes two one-dimensional harmonic oscillators with a cubic inter- 
action. It is one of the simplest Hamiltonians to display soft chaos in classical mechanics: by increasing 
the total energy a transition from an integrable to an ergodic system is induced. Originally conceived to 
model the chaotic motion of stars in a galaxy it later became an important milestone in the development 
of the theory of chaos [|3|], partly because of the conceptual simplicity of the model. 

In order to investigate this continuous loss of integrability with growing total energy Gustavson 
transformed the classical Hamiltonian by a series of canonical transformations to a Birkhoff-Gustavson 
normal form [^], which allowed him to construct an additional constant of motion. Thus he was able to 
analytically reproduce Poincare surfaces of section as obtained by numerical integration. 

The present work is based upon two publications ([^ and [Q]) in which Delos and Swimm used 
the classical Birkhoff transformation [||] as applied by Gustavson to analyze how the classically chaotic 
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behavior of the system is transformed into quantum mechanics. The Birkhoff-Gustavson normal form is 
a power series in oscillator Hamiltonians and thus allows a direct determination of a quantum mechanical 
eigenvalue spectrum from the classical Hamiltonian. For a fixed set of parameters Delos and Swimm 
analytically calculated the spectrum of the Henon-Heiles Hamiltonian which reproduced the eigenvalues 
obtained by numerical diagonalization of finite matrices. 

The main problem of the quantum mechanical Henon-Heiles Hamiltonian is the fact that it is not 
bounded from below. It contains a cubic potential. The classical motion discussed by Gustavson and 
by Delos and Swimm corresponds to initial conditions near the local minimum of the potential and to 
an energy that is below the saddle point value of the potential. By these conditions the classical motion 
is always restricted to a finite region. In the corresponding quantum problem, the particle will always 
tunnel through the barrier. Therefore the eigenvalue spectrum calculated by Delos and Swimm is not the 
real eigenvalue spectrum of the Hamiltonian. It describes effective states that can be used to describe the 
dynamics near the minimum of the potential and for times that are small compared to the escape time. 

The aim of the present paper is to first quantize the classical Hamiltonian and to then diagonalize 
it using the method of flow equations which was introduced by Wegner [Q] in 1994. It is clear that 
concerning the tunneling problem, the flow equations have the same limitation as the quantization of the 
Birkhoff-Gustavson normal form by Delos and Swimm. The bound states and the eigenvalues obtained 
using flow equations allow only an effective description for small times. One advantage of the flow 
equations is that a correct and simple treatment of the system is possible even if the two frequencies of 
the harmonic oscillators are commensurate. This is not the case if Birkhoff-Gustavson normal form is 
quantized. Furthermore, the quantum mechanical treatment allows a description of the dynamics. 

The structure of the paper is as follows. The next section offers a general introduction to the flow 
equation method. In sections ^ and ^ their application in two variants to the Henon-Heiles model is 
treated. Section § contains some results of the diagonalization: in a table the energy eigenvalues obtained 
in the flow equation approach are compared with those obtained by numerical matrix diagonalization for 
a fixed coupling constant. A graph shows the dependence of the calculated eigenvalues upon the coupling 
strength. In section |6| a case of commensurate frequencies is treated and a similar table of eigenvalues 
is obtained. Section ^ offers a method to investigate the dynamics of the quantum mechanical system, 
transition amplitudes between the eigenstates of the uncoupled system are determined. The effect of 
growing coupling strength upon the transition amplitudes is shown. The last section contains a summary 
of the results of the present work, gives a comparison to the work of Delos and Swimm and a discussion 
of the limitations of the Birkhoff-Gustavson transformation. 



2 Flow Equations 

The method of flow equations consists in a continuous unitary transformation of a given Hamilton oper- 
ator H, which can be written in differential form: 

^ = [Ti(^),//(£)]. (2.1) 

There are several possibilities to choose the antihermitean generator r\ so that //(oo) becomes diagonal. 
Wegner [^] proposed: 

ri(^) = [Hd{t),H{l)] = [Hd{l),Hr{t)] , (2.2) 

where and H, are the diagonal and the off-diagonal portions of the Hamilton operator respectively. 
A detailed argument for the usefulness of this choice of the generator can be found in [Q]. But the 
consistency can be easily verified, since in the limit £ ^ oo as //(£) becomes more diagonal r[{t) will 
vanish and so will ''^jp . 
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The flow equation method has been apphed to various models (see e.g. [|, |, |I0[ |ri| |12|]). The general 
behavior is such that terms that appear in r\ < ^J^ will result in new terms in the transformed Hamilton 



operator by (|2JJ). If this iterative process does not result in a closed set of differential equations it can be 
forced into such by defining an order for the appearing terms - e.g. the number of creation operators in 
them - and neglecting all terms of higher order. This approach will be called cut-off. 



A second approach to handle the system of differential equations (2.1) and (2.2) for a given Hamilton 
operator H = //^(O) +'kHr{0) is to define the transformed Hamilton operator as a power series in the 
coupling constant X 



(2.3) 



k=Q 



The coefficients Hk{i) can be determined iteratively and this method shall be called iteration. Both of 
these procedures - the cut-off and the iteration - will become more transparent as they are applied to the 
Henon-Heiles Hamiltonian in the next two sections. 




Figure 1 : Henon-Heiles potential ( |3.lD 
(X = -0.1, ?i = 0.1, w = 1.3, v = 0.7) 



3 The Cut-off Procedure 

The Henon-Heiles Hamiltonian can be expressed as a function of two spatial coordinates qi and q2 and 
the two momenta p\ and p2 (see figure [l]): 

H = ^-w{p\ +q\) + K{pI + ql) + Xq2{ql + nql) . (3.1) 
Quantizing this classical Hamiltonian using the operators 

a:=^{q + ip), a' := ^{q - ip) . (3.2) 
will give the Hamilton operator 

H = wa^a + vb^b + l {b^ + b) {{a^ +af + n [b^ + bf) , (3.3) 
where the coupling constant X has been rescaled and the constant term dropped. 



This Hamilton operator is to be transformed into tlie quantum mechanical equivalent of a Birkhoff 
normalform, a power series in oscillator Hamiltonians: 



k,m=0 



(3.4) 



k.m=0 



The generalized frequencies wtm are obtained from the transformed Hamilton operator H{£) in the limit 
/ ^ oo. A consistent ansatz is H{£) = Hd(i) +Hr{tj where: 



woo(^) + W2o{t) a^^a^ + W02 {£) b^^b^ + wn{i) abh 



and 



(0^2 + (^t + ^) + (^12 _ ^2) (^t _ ^) ^2 (i) 

+ ia''a){b'f+b)x^{e) + (Z7l-3+Z73)x4(£) 
+ {b^^b + b^b^)xs{l) + (Z?^+ft)x6(^) 



The ^-dependent coefficients are determined from (2.1) and ([2.2|) which combine to: 



(3.5) 



(3.6) 



Neglecting all terms of third order in X or higher we obtain the following set of differential equations by 
coefficient matching: 



w'oo = 


-4x?v - 


Sx^w — 8X1X2V — \6x1X2w — 36x4V — 


9 9 9 
2X6V — 4X2V — 8X2^ 


w' = 


-8x?v - 


ISxfw — 1 6x1x2 V — 32x1X2^ — 2x3 V - 


9 9 

- 4X3X6V — 8x2V — \6X2W 


1 

V = 


— \6x\w 


— \6x\X2v — IO8X4V — 4X5V — 8X5X6V 


— \6X2W 


W2Q = 


-Ax\v - 


\6x\X2W — 2x3 V — 4X2V 




w'n = 


—32x\w 


— 32x\X2V — 8X3X5V — 32X2H' 




w'o2 = 


-54x4V 


-6x^v 




x\ = 


— XlV^ — 


4x\w^ —4x2vw 






—Ax\vw 


— X2V^ — 4X2 






—x^v^ 






f 

— 


-9X4 v^ 








-X5V^ 






/ 

^6 = 


-X6V^ 







(3.7) 

The equations for X3, . . . ,X6 are all of similar type, the equations for w,^ are uncoupled. The system can 
therefore be reduced to five equations. It could not be solved analytically. But the asymptotic behavior 
for large I is the following: 

Assume that w{l) « w{°°) = w<x> and v{t) v{°°) = Voo. Then the off-diagonal elements show an 
exponential decay (except in the case of commensurate frequencies 2woo + Voo = 0): 

x\{tj = coexp (-(2w«> + Voo)^ 
X2{tj = coexp (-(2wt>o + Voo)^ 



(ioexp I 



(3.8) 
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The diagonalized Hamilton operator was determined by numerical integration of ( ^?7| ) using a Runge- 
Kutta procedure. For various sets of fixed parameters the approximation was improved by extending the 
calculation to all terms up to fourth order in X. The resulting set of differential equations for coefficients 
xi , . . . ,X48,w, V and the coefficients Wij corresponding to the diagonal operators {a^'a') {b^^y) were also 
determined by numerical integration. An eigenvalue spectrum obtained from the transformed Hamilto- 
nian in the limit C ^ oo was calculated for different values of the coupling constant X. Results of these 
calculations are presented in section ^ 



4 The Iterative Procedure 



A far more elegant way to solve the flow equation (gjj) for the Henon-Heiles Hamiltonian is an iterative 
calculation of the Hk{(-) defined in (2.3). It avoids the numerical integration applied in the previous 
section and allows a better insight into the transformation mechanism, specifically the behavior in the 
case of commensurate frequencies. 

The Hamiltonian ( |3.3| ) is given in the form 

H{<d)=HQ{Q)+XHi{Q) with Ho{Q) =wa^ a + vb'^b . (4.1) 



The transformed Hamiltonian is defined as a power series in X (23). In deviation from the original choice 
(p!2[) the generator is now defined as: 



^{l):=[Ho,H{t)] = Y,X\k{e), where ^k = [HM. (4.2) 

This choice of r\ is in accordance with the classical Birkhoff transformation as applied by Gustavson 
It makes the iterative calculation simpler, since commutation with Hq will reproduce a given operator 
ttrnva^^db^'^b": 



Ho, a' a 



etrmna'^'a'W, (4.3) 



where Ekrmn ■= [(^ — + (m — n)v] . Inserting the power series (23) for the transformed Hamiltonian 
into the flow equation ( |2.l[ ) and comparing the coefficients of the powers of X, gives the differential 
equations: 

^^^M = [[Ho,H„],Ho] + [[Ho,Ha],Hh] V« = 0,l,2.... (4.4) 
It follows that Hq = const, and Hk{£), k = 1,2 . . . can be iteratively calculated: inserting the general ansatz 

H„{e)= ^ 5krmn{i)a'^W'%" (4.5) 

k,r,m,n 

into ( |4.4| ) results in differential equations of the form 

-J^hnnn{i) = -^krmnhrmn{t) + Cf-kmini^) (nO Summation) , (4.6) 

where Zkrmn > was defined above. The function (Xurmn (^) can be shown to be a sum of terms of the form 
crexp(-Y^) (c = const., Y>0, ?i = 0, 1,2. . . ) . (4.7) 
The solution to ( |4~6| ) is 

hrmni^) = SXpi-Ekrmn ■ d l' akrmn{f ) ^'^Vi^krmJ' ) ■ (4.8) 

JO 

Ignoring for a moment the inhomogeneity CLurmni^) we find the following behavior: all terms will decay 
exponentially except for two cases in which Zkrmn may vanish: 
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1. The given operator term a^'^a''b^'"b" is diagonal (k = r and m = n). 

2. The frequencies w and v are commensurate. 

It can be seen that the inhomogeneity CUkrmni^) will not change this general behavior. Thus the trans- 
formation can be successfully performed in the case of incommensurate frequencies and only diagonal 
terms will remain in the limit £ — > oo. In the case of commensurate frequencies off-diagonal operators 
a^'^a''b^'"b" may only remain for Ekrmn = 0. 

Since the number of terms in Hk{£) grows rapidly with k the described procedure with all its alge- 
braic manipulations was performed by a computer program in the language C. This way Hq,... were 
determined. The results of these calculations will be presented in the following section. 



5 Results of the Diagonalization 



The two procedures described in the last two sections were applied to the Henon-Heiles Hamiltonian ( |3.1| ) 
in a case of incommensurate frequencies with the values w = 1.3, v = 0.7, X= —0.1, n = 0.l. The iteration 
procedure to 8th order results in the following diagonal form (at X = —0. 1): 

H = 2.20910 • 10"'a'''^£?^ + l. 19855 •10"'^a^'*aVZ; + 2. 17973 •lO^^a^'^ a'* 
+ 1.54236- lO^^a^^ a^Z^^^Z^^ +4.12440- 10" a^b^b 
-8.65783 - IQ-^a^^ + 6.09731 - 10'^ a^^ a^b^^b^ 
-4.99568 - lO^^a^^ a^b'^^b^ - 0.00030a'^^ a^b^b - 0.00634a^2a^ 
-2.46444 - lO^^a^ ab'^'^b'^ - 6.97206 - lO^^a^ ab'^^b^ - 0. 00022a^ ab^^b^ 
-0.01 Ula^ab'^b+ 1.28264a^a + 3.44234 - lO'^b^^ b^ 
-2.84772 - 10"^ b^^b"^ - 1.59258 - lO'^b'^^ b^ -0.0011 Ib'^^b'^ 

+0.69148Z^^Z^ + 0.99552 (5.1) 

Table [l| shows the eigenvalue spectra derived from the transformed Hamiltonian for the cut-off pro- 
cedure carried out to fourth order and for the iteration procedure to fourth and to eighth order in X. The 
numerical calculations were performed by diagonalizing a 900x900-matrix using FORTRAN-routines 



presented in []130. The error A gives the relation of the deviation from the numerical value with respect to 

(5.2) 



the total shift due to the coupling: 



''num. ^frc: 

It was calculated for the values of the iteration procedure to 8th order. 

The calculated spectrum reproduces the numerical data well. For the fixed value of X = 0.1 the 
iteration procedure gives better results than the cut-off procedure. Iteration to higher orders does not 
necessarily improve the approximation of the numerical data: the ground state energy is more accurate 
for the iteration to 4th order than for the iteration to 8th order. Therefore there is no monotonous con- 
vergence of the determined eigenvalues to the exact ones with increasing order of iteration. Whether 
the eigenvalues determined in the flow equation procedure are always above the exact ones could not be 
shown analytically. 

One can determine how the energy eigenvalues for a given pair of quantum numbers changes as 
the coupling strength X is increased. The result is shown in figure ^ for 12 eigenvalues in the cut-off 
procedure to 3rd order and the iteration procedure to 8th order. With growing coupling X the potential 
well in figure [T] becomes shallower and the eigenstates move closer together. The eigenvalues decrease as 
the coupling is increased. Moreover, one finds that states with higher energy at X = —0.1 will drop faster 
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Table 1 : Comparison between flow equation calculations and numerical data in a case of 
incommensurate frequencies (w= 1.3,v = 0.7,X= — 0.1,« = 0.1) 



n 




«2 


Cut-off 


Iterat. 4 


Iterat. 8 


Numerical 


^free 


A(%) 


1 


n 
u 



u 


V.yyJJO 1 


QQ'^'^9 1 


U.yy J JZJ 


yj.yyJJ Ly 


1 00 


1 ^^8QQ 


9 



u 


1 

i 


1 (^87949 
l.Do 1 L^L 


l.Oo 


1 687010 
l.Do /UlU 




1 70 
1. /U 


1 9^090 


J 


1 

1 


A 

u 


9 97S'^4'? 
Z.Z / OJ^O 


9 9781 70 
Z.Z lolly 


9 978 1 70 
Z.Z / o 1 /U 


9 9781 ^^9 
Z.Z / o 1 jZ 


9 ^O 
Z. jU 


O 1 7^770 

yj.i 1 J 1 lyj 


4 





2 


2.375702 


2.375106 


2.375064 


2.375036 


2.40 


0.112162 


5 


1 


1 


2.959696 


2.958536 


2.958439 


2.958353 


3.00 


0.206497 


6 





3 


3.060918 


3.059734 


3.059592 


3.059551 


3.10 


0.101362 


7 


2 





3.549267 


3.548183 


3.548119 


3.547947 


3.60 


0.330432 


8 


1 


2 


3.637534 


3.635141 


3.634827 


3.634664 


3.70 


0.249480 


9 





4 


3.742857 


3.740832 


3.740491 


3.740435 


3.80 


0.094015 


10 


2 


1 


4.219694 


4.216860 


4.216555 


4.216180 


4.30 


0.447387 


11 


1 


3 


4.312030 


A. 301916 


4.307197 


4.306912 


4.40 


0.306162 


12 





5 


4.421492 


4.418334 


4.417653 


4.417578 


4.50 


0.090995 


80 


1 


14 


11.502109 


11.437765 


11.386826 


11.348431 


12.10 


5.108646 


81 


8 


1 


11.496693 


11.460618 


11.465128 


11.412886 


12.10 


7.603105 


82 


7 


3 


11.535246 


11.478808 


11.477820 


11.415802 


12.20 


7.908462 


83 


6 


5 


11.591118 


11.517846 


11.506570 


11.432484 


12.30 


8.540015 


84 


5 


7 


11.664419 


11.578700 


11.552905 


11.470273 


12.40 


8.887770 


85 





16 


11.659609 


11.614802 


11.583460 


11.532429 


12.20 


7.644281 



than the lower states as the coupling is increased. This is due to the fact that the eigenstates to higher 
eigenvalues are more spread out in space such that the effect of the X • ^j-term in the potential {3A) upon 
them is larger. 

Comparison with numerical data seems to indicate that for larger coupling values of X the cut-off 
procedure gives more accurate results than the iteration procedure. However, a precise quantitative anal- 
ysis is not possible in the X-range in which the eigenvalues from the two flow equation procedures differ: 
for growing values of X one has to restrict the numerical diagonalization to smaller matrices in order to 
avoid the effect of the continuum causing the appearance of intermittent states (seen in figure ||). 

6 A Case of Commensurate Frequencies 



The two flow equation procedures (sections ^ and ^ were applied to the Hamiltonian ( 3J_ ) in a case of 



commensurate frequencies (w = 1.0, v= 1.0,?i = — 0.1,w = 0.1). As shown in section off-diagonal terms 
will not entirely disappear from the transformed Hamiltonian due to the commensurability. According to 
( p3| ) and ( |4.6D remaining off-diagonal operator terms are of the form a^'^a''b^'"b" with Ekrmn = {k — r) + 
{m — n)= 0. They couple states | «i , «2) for which «i + ?i2 = const. Neglecting these off-diagonal terms, 
we obtained fairly accurate eigenvalues. To account for the off-diagonal terms small tridiagonal matrices 
were numerically diagonalized within the originally degenerate subspace. This improves the precision 
of the calculated eigenvalues. The cut-off procedure does not have this problem. Since the oscillator 



frequencies v and w in (3.3) are £ dependent (see (3.7)), initially commensurate frequencies become 



incommensurate for finite i. The results for the cut-off procedure to fourth order, iteration procedures to 



7 



5.0 
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3.0 



2.0 
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0.0 




, \ , \ 

-0.34 -0.24 -0.14 



X 

Figure 2: The eigenvalues for the quantum numbers of the lowest 12 
eigenstates at X = —0.1 as a function of the coupling strength X in the 
Cut-off procedure to 3rd order (above) and the iteration procedure to 8th 
order (below). The numerical integration breaks down for large coupling 
(above). 



4th and 6th order and the improved values of the 6th-order iteration are listed in table g. A gives the error 
of the improved values with respect to the total shift due to the coupling: 



-'impr. 



^num. Efree 

The originally degenerate subspaces are separated by horizontal lines. 



(6.1) 



7 Dynamics of the Quantum Mechanical System 

In the framework of classical mechanics one can put a particle in the potential well (figure [T]) and calculate 
the trajectory for a given set of starting values. To model this classical approach in the framework of 
quantum mechanics one can ask: How does a state originally located within the potential well evolve 
with time ? The eigenstates of the uncoupled oscillator are located within the potential well. Their time 
evolution is given by the Hamilton operator H of the coupled system. The absolute value of the matrix 
element 

(p|exp(///0l«) (7.1) 
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5.5 



4.5 



3.5 



2.5 



1 .5 



A. 

Figure 3: Numerical diagonalization of a 900x900-matrix. For 256 
X-values in the interval [-0.5,-0.1] the first 29 eigenvalues E above 0.7 
were determined. With growing perturbation the effect of the continuum 
becomes dominant. 



indicates how much of a particle is in state | P) after time t if the particle was located in state | a) at time 
where | a) and | (3) represent two of the eigenstates of the uncoupled oscillator. 

To calculate such matrix elements the eigenstates of the uncoupled oscillator are expressed in terms 
of the eigenstates of the full Hamiltonian of the coupled system. For this purpose one can set up flow 
equations for the transformation of states which results in a fairly large set of differential equations. 
Since the transformation of the Hamiltonian was already calculated there is an easier way to calculate the 
transformation of states which will be explained in the framework of the iteration procedure introduced 
in section 0: 

In analogy to the transformation of the Hamiltonian one determines the transformed annihilation 
operator a{€) = U\£)aU{£) from the flow equations 

^ = m,m, (V.2) 

and the generator r\ already known from the calculation of H{£). One obtains the annihilation operator 
as a power series in the coupling constant X similar to the one obtained for the transformed Hamiltonian 
p^ ). By definition the transformed ground state of the uncoupled oscillator |0) is given by: 

a(£ = oo)|0)=0 , b{£ = oo)\0)=0 and (0|0) = 1. (7.3) 

These equations can be solved for |0) as a function of the eigenstates \n,m)^ of the full Hamiltonian: 

oo 

|0)= £ SnX'\n,m)^. (7.4) 

n,m,i=0 

One can then construct any exited state | a) = a^^ (oo) b'^'" (oo) 1 0) . In the following the abbreviations 
a: = a{oo) and b:=b{oo) will be used. 
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Table 2: Comparison between several flow equation calculations and numerical data in a 
case of commensurate frequencies 

iw = \.0,v=l.O,X = -0.l,n = 0.\) 



n 


«i 


«2 


Cut-off 


Iterat. 4 


Iterat. 6 


improved 


Numerical 


Efree 


A(%) 


1 








0.997021 


0.996989 


0.996990 


0.996990 


0.996987 


1 


0.084 


2 


1 





1 983815 


1.983444 


1 983415 


1 983415 


1 983420 


2 


0.030 







1 


1 991400 


1 991 1 87 

J- . W J- J- O / 


1 991 1 80 

J- . W J- J- o\j 


1 991 1 80 


1 991 170 

J- . W J- J- 1 \J 


2 


0.11 1 


4 


2 







2 961007 


2 9608 1 6 




2 957081 




195 




1 

1 


1 

1 


9 968458 


2 9671 1? 


9 966985 


2 966985 


2 966957 

Z.. y\j\jyo 1 


'X 


084 


6 





2 


2.985028 


2.984515 


2.984469 


2.988288 


2.988269 


3 


0.162 


7 







3 932391 


3 929514 


3 928861 

VZ.OOVJ J. 


3 917954 


3 91 8277 

O . y 1. oz* / / 


4 


395 


8 


2 


1 


3.937177 


3.933515 


3.933009 


3.926540 


3.926527 


4 


0.018 


9 


1 


2 


3.952348 


3.949675 


3.949312 


3.960219 


3.960177 


4 


0.105 


10 





3 


3.977903 


3.976969 


3.976841 


3.983311 


3.983255 


4 


0.331 


11 


4 





4.894173 


4.888802 


4.887144 


4.865409 


4.865039 


5 


0.274 


12 


3 


1 


4.897560 


4.890232 


4.888878 


4.870960 


4.871053 


5 


0.072 


13 


2 


2 


4.911331 


4.904683 


4.903555 




4.916403 


5 


2.352 


14 


1 


3 


4.935486 


4.931130 


4.930349 


4.948267 


4.948041 


5 


0.435 


15 





4 


4.970025 


4.968546 


4.968277 


4.979130 


4.978266 


5 


3.973 


16 


4 


1 


5.849605 


5.837100 


5.834141 


5.796496 


5.795370 


6 


0.550 


17 


5 





5.847617 


5.838708 


5.835184 


5.799784 


5.799166 


6 


0.308 



For several final states | p ) matrix elements of the type (^) are determined as a function time. Be- 
cause of the coupling of time and energy in the exponential function an expansion of the exponent in 
powers of X - and thus in powers of t - does not give the right long term behavior. Therefore the expo- 
nential functions are not expanded. The resulting matrix elements are of the form: 

{^\exp{iHt)\a)=Y,at^'^e'^^', (7.5) 

k 

with coefficients and exponents b^. represents a difference of energy eigenvalues which are deter- 
mined in the iterative process to the 6th power of X. Thus if the coefficients of the expansion in X are all 
of order 1 the phases E^t can be determined up to 1 percent accuracy in the range 

t < t,nax ~ ^ = 10^ (7.6) 

This gives only a rough estimate since Siegel [ |T4{ ] proved that the Birkhoff normal form will generally 
not converge. This can easily be deduced from the fact that the Henon-Heiles potential is not integrable 
whereas any polynomial in Birkhoff normal form will always be integrable. One expects an asymptotic 
convergence of the expansion in X such that the coefficients of higher powers of X will generally grow. 
Thus the time range in which the calculated matrix elements are valid is smaller than the one given above. 

For the initial state | a) = a ' 1 0) amplitudes for the transition to the six final states | P) = 1 0) , a^Z^ ' 1 0) , 
a^b^^\0), a^^\0), a^^b^O) and a^b^^\0) were determined at X = —0.1. Figure ^ shows the first four of 
them as a function of time. Amplitudes for transitions to final states with an even number of a^- operators 
can be shown to vanish. 
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Figure 4: Transition amplitudes f{t) = {0\aUa^\0), {0\abUa^\0), 
{0\ab'^Ua^\0) and {0\a^Ua'^\0) (in order of size) where 
U := exp{iHt), a = a{oo) and b = b{°°) 



The sum of their absolute values was subtracted from 1 - in figure § - to show that these are indeed 
the relevant transitions. The amount of negative value in this figure gives an indication of the numerical 
error. 



The parameter X is a measure of the strength of coupling between the two harmonic oscillators (3.1). 
One expects that a growing coupling value facilitates transitions between different states. Figure |6| shows 
the square of the transition amplitude (0|aexp(///f)a^|0) as a function of time for various values of X. 
Moreover, one finds that the dominant oscillation frequency decreases with growing coupling. 

In the framework of quantum mechanics one expects a particle located within the potential well of 
figure [T] to tunnel through the potential barrier with a certain probability. This probability depends upon 
the energy of the particle and upon the size and geometry of the barrier. Thus the tunneling should 
become apparent if one investigates the dynamics of higher exited states or if one decreases the size of 
the potential well by increasing the coupling X. Estimates of the tunneling frequency show that the effect 
of tunneling is negligible for X = —0.1. However, for X = —0.25 tunneling should become apparent. In 
figure ^ one finds a modulation in the amplitude of oscillation for X = —0.25. But this may also be a 
numerical artifact because in this range of coupling strength the numerical precision decreases. Since 
the Birkhoff approach consists in approximating a power series by a finite polynomial, any transition 
amplitude will always be a finite sum of harmonic oscillations. Therefore any calculated time evolution 
will necessarily be periodic in time. 

Physical observables can be calculated in essentially the same way as was shown for the transition 



amplitudes (7.1). Using the relations 



1 

7! 



(a' +a) 



V2 



(7.7) 
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0.0010- 




FigureS: (l -LpKPI^^I«)P) with initial state |a)=al'|0) and 
the six final states | P) mentioned above (U := exp(///f)) 




0.0 J ' < . < ' < ' < . 1 

0.0 10.0 20.0 ^ 30.0 40.0 50.0 



Figure 6: Evolution of the first exited state 

for the coupling values ?i= —0.1, —0.15, —0.2 and —0.25 

(in the order of growing amplitude), where f{t) = 

{0\aexp{iHt)a^\0). 



for spatial and momentum coordinates one can easily calculate expectation values of the form (q^) = 
{a\U^q^U\a) , which show a similar oscillatory behavior as that in figure ^. 
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8 Summary and Outlook 



In this paper we calculated effective eigenvalues of the quantum mechanical Henon-Heiles Hamiltonian 
using flow equations - a method of continuous unitary transformation proposed by Wegner We used 
two different procedures to solve the flow equations - an iterative procedure and a Cut-off procedure. The 
Cut-off procedure has been used before in several applications of the flow equations. It seems to be most 
appropriate if a (perturbative) renormalization of the Hamiltonian has to be done. This is not necessary 
in the present case. The Cut-off procedure has the disadvantage that in most cases it is difficult to solve 
the resulting differential equations explicitly. Often one can extract the asymptotic behavior and based 
thereon an approximate solution. But if one wants to have precise numerical results, one has to solve the 
differential equations numerically. The main advantage of the iterative procedure is that the differential 
equations can be solved explicitly and that it can be carried out to much higher orders. 

We used both methods to calculate the eigenvalue spectrum of the Henon-Heiles Hamiltonian. In 
a case of incommensurate frequencies the eigenvalues coincide well with those obtained by numerical 
diagonalization of a finite matrix. In the treated case the precision of results is comparable to that obtained 
by Delos and Swimm who used consecutive canonical transformations to approximate the classical 
Hamilton function by a Birkhoff normal form. For small coupling value the iteration procedure gives more 
accurate results whereas for larger coupling values the cut-off procedure seems to be better. However, 
for very large coupling the potential well decreases in size and the discrete spectrum disappears. This is 
found in a series of numerical matrix diagonalizations at various coupling values. 

The case of commensurate frequencies can be treated in essentially the same way, matrix diagonal- 
ization shows good agreement with numerical results. In the commensurate case the diagonalization 
procedure is much less lengthy than that presented by Delos and Swimm. A quantitative comparison to 
their results is not possible since the normal form cited in does not correspond to the given parameter 
set. And the presented eigenvalue spectrum does not match with either the parameter set given or the 
normal form cited. 

On the basis of the transformation of eigenstates we analyzed the quantum mechanical dynamics of 
the Henon-Heiles system. We determined the time evolution of a set of states located within the potential 
well. We found oscillations between different states with amplitudes depending on the coupling strength. 
A certain completeness is found by adding up absolute values of several transition amplitudes. 

It is interesting to see that continuous unitary transformations can be used to obtain precise results 
for a system like the quantum mechanical Henon-Heiles model. The method was originally designed to 
calculate eigenvalues of a given Hamiltonian, at least approximatively and close to the ground state ||^. 
In many applications of the method, the goal was to obtain an effective Hamiltonian that describes well 
the low energy behavior of the system (see e.g. [^, ^ [I^, [ri|]). The effective Hamiltonian calculated in 
the present approach is in a quantum mechanical Birkhoff normal form, i.e. a polynomial in the oscillator 
Hamiltonians {p^+q^). It does not describe the ^^-shape of the Henon-Heiles potential for large absolute 
values of q. The energies we calculated using this method are energies of approximately stationary states. 
As a consequence one cannot describe the tunneling through the barrier. This limitation is not imposed 
by the flow equation method. Even numerical diagonalizations will not allow a correct description of 
the tunneling. But for sufficiently large escape times the results describe the short time behavior of the 
system quite well. 

The transformation to a Birkhoff normal form yields an asymptotic series because the classical 
Henon-Heiles system shows a transition to a chaotic regime. In principle the flow equations can be 
successfully applied to a quantum mechanical Hamiltonian that has a chaotic classical counterpart (see 
also [12]). A first focus of interest for such a system are statistical properties of the spectrum [^. Flow 
equations can be used to calculate precise eigenenergies. The way normal ordering is introduced deter- 
mines the interval of energy where precision is high. Typically one has to restrict the number of couplings 
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in the Hamiltonian using some truncation scheme. The neglected terms are in a normal ordered form. If 

one chooses the ground state as the basis for the normal ordering, the resulting effective Hamiltonian 
describes the ground state and low lying excitations quite well. In the present approach the normal order- 
ing was introduced with respect to the ground state of the uncoupled system. Therefore we were able to 
obtain precise results for energies close to the minimum of the potential. One could also choose a normal 
ordering with respect to some high energy state in order to calculate eigenenergies close to this energy 
scale. While this approach is not sensible in the case of the Henon-Heiles system, because the potential 
well has only a limited depth, it can be useful for other physical systems. 
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